Competing superconducting, magnetic and charge orderings in 
the AF Heisenberg-Kondo lattice with Dirac electrons 
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Abstract. - Many recently discovered advanced materials, such as high-Tc cupratcs, iron pnic- 
tides and several heavy-fermions, exhibit a rich phase diagram suggesting the presence of different 
competing interactions that would lead to various types of ordering. Nevertheless, there is not yet 
a clear unifying picture allowing the understanding of the detailed mechanisms that generate such 
competing interactions. Blaving such a picture, however, could quite well be at the very roots of 
the requirements for understanding high-Tc superconductivity in cuprates and pnictides, for in- 
stance. In this work we consider the antiferromagnetic (AF) Heisenberg-Kondo lattice, consisting 
of localized spins with AF exchange interactions between nearest neighbors on a square lattice 
and itinerant electrons, which undergo a magnetic Kondo interaction with the localized spins, 
but are otherwise non-interacting. Using the Schwinger-boson (CP^) formalism and assuming the 
electrons are Dirac-like, we integrate on the localized degrees of freedom thereby obtaining the ef- 
fective interaction among the itinerant electrons. This contains a BCS-like superconducting term, 
a Nambu-Jona-Lasinio-like, charge gap term and a Ising and Heisenberg-like magnetic terms. 
All these four competing interactions, therefore are generated by the original Kondo magnetic 
interaction. 



Introduction. — Many recently discovered systems in 
' condensed matter physics present a deep interplay among 
, different types of orderings such as superconducting, mag- 
' netic or charge ordering. Also in some of them the elec- 
tronic excitations behave as Dirac fermions These include 
iron pnictides, superconducting cuprates, transition metal 
dichalcogenides (TMD) and possibly graphene, if the ma- 
terial is doped, strained or has atoms adsorbed on its sur- 
face. 

Iron based pnictide materials, for instance, undergo a 
transition from a magnetically ordered state to a super- 
conducting one upon doping [IH3] . For the particular case 
of the 122 materials, magnetic order and superconductiv- 
ity coexist in a small region of the phase diagram j4] and 
the new quasiparticles in that region exhibit a Dirac-like 
linear energy dispersion [SHZ] with Cooper pairing possi- 
bly promoted by short wavelength antiferromagnetic spin 



fluctuations. 

On the same token, the parent compounds of the 
cuprate superconductors are insulators presenting AF or- 
der. As charge carriers are added to the Cu02 planes, 
there is the onset of superconductivity, with the character- 
istic dome-shaped superconducting phase diagram [TOHTT] . 
Dirac points appear in the intersection of the nodes of 
the d-wave superconducting gap and the two-dimensional 
(2D) Fermi surface |T8H20) . Strongly interacting 2D Dirac 
fermion systems also exhibit a dome structure in their su- 
perconducting phase diagram [HI [35], so we may spec- 
culate whether Dirac fermions may play any role in the 
description of the cuprate superconductors. 

Moroever, the quasi 2D TMD are layered compounds 
where s-wave superconductivity coexists with a charge 
density wave (CDW) at low temperatures and applied 
pressure [5HlO). A theory has been proposed in which 
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Dirac electrons appear close to the nodes of the CDW gap 
and form Cooper pairs [TT1[T2]. The theory is consistent 
with the linear decay of the temperature dependent critical 
field [14] ■ which is observed experimentally in the copper 
doped dichalcogenide Cua;TiSe2 [TS]. 

Graphene, on the other hand, is a semi-metal with gap- 
less electronic quasiparticles, which due to the peculiar 
lattice structure behave as Dirac fermions. Pure graphene 
does not exhibit superconductivity or magnetism, how- 
ever, upon doping and/or straining may display competing 
orders, such as local magnetic moments, superconductiv- 
ity or an excitonic gap |23| . 

Heavy fermions is another vast class of materials pre- 
senting a rich competition of different types of order in 
their phase diagrams [24H26] 

The above mentioned systems have been exten- 
sively studied both theoretically and experimentally but, 
nonetheless, we still do not have a clear unified picture af- 
fording a detailed understanding of the microscopic mech- 
anisms that lead to each kind of ordering. It would be very 
instructive and inspiring, therefore, to find fully control- 
lable models for systems displaying phases with the afore- 
mentioned different types of ordering, where one could pre- 
cisely trace back the original interaction and the mecha- 
nisms leading to such phases. 

In this work, we investigate the AF Heisenberg-Kondo 
lattice model [27| , a system containing both localized spins 
and itinerant electrons with a mutual Kondo-like magnetic 
interaction. Our main aim is to determine what is the 
effective net interaction among the conduction electrons, 
which results from their magnetic interaction with the AF 
substrate of localized spins. For this purpose, we use the 
well-known Nonlinear Sigma Model (NLSM) description 
of the latter [35]. The itinerant electrons, conversely, are 
supposed to have a tight-binding band structure showing 
the presence of Dirac cones whose vertices touch at the 
interface between the valence and conduction bands. We 
assume the system to be close to half-filling and there- 
fore describe the kinematics of the itinerant electrons by 
the Dirac hamiltonian. In order to describe the the mag- 
netic Kondo interaction it is convenient to re-phrase the 
NLSM in the CP^ language [29], whereby that interaction 
becomes a quartic term involving two Dirac fermion fields 
and two bosonic (Schwinger Boson) fields. 

Our strategy will be to functionally integrate over the 
bosonic fields in order to derive the resulting effective in- 
teraction existing among the fermion fields. In this pro- 
cess, we show that the Kondo interaction among itiner- 
ant and localized electrons can be completely expressed 
as a gauge coupling between the Dirac fermions and the 
Schwinger bosons, mediated by the CP^ vector gauge field. 
Our final result for the effective electron interaction is then 
obtained upon integration over this gauge field. 

The resulting interaction possesses four pieces: a super- 
conducting BCS-type term, an Ising and a Heisenberg-like 
magnetic interactions and a Nambu-Jona-Lasinio-type in- 
teraction. These terms will favor respectively supercon- 



ducting, magnetic and insulating charge-gapped ordering. 
We conclude that the original system of localized and itin- 
erant spins yields ultimately an interacting electronic sys- 
tem with those types of competing orders. 

The Model and Its Continuum Limit. We con- 
sider a single layered system containing both localized and 
itinerant spins in which the former are located at the sites 
of a square lattice and have an antiferromagnetic exchange 
integral while the latter are conduction electrons with a 
tight-binding dispersion relation, which is assumed to be 
Dirac-like. The localized spins mutual interaction will 
be described by an AF Heisenberg hamiltonian on a 2D 
square lattice, whereas their interaction with the itinerant 
ones, by a Kondo-like term. The complete hamiltonian, 
therefore, contains three parts, namely 



H = J^Si-Sj-tY^ (cLcja + h. 



(iJ) 



(ij) 



(1) 



where S; is the localized spin operator and cj^ is the cre- 
ation operator of an itinerant electron of spin a =t j 
both at site i. Frequently we have materials for which 
there are electrons coming from different bands or even 
from inequivalent regions of the Brillouin zone. In these 
cases we would add an extra label a = 1, to the elec- 
tron operators. Here, for the sake of simplicity, we shall 
omit such a label, this fact having no effect in our conclu- 
sions. 

In order to obtain the partition function, we employ the 
continuum path integral approach. By using a basis of spin 
coherent states we have the localized spins Si replaced by 
their correspondidng eigenvalues: S'N(x), where 5* is the 
spin quantum number and |N(x)p = 1. N is then decom- 
posed into two perpendicular components associated re- 
spectively with ferromagnetic and antiferromagnetic fluc- 
tuations 



N(x) = aL(x) + (-l)l''Vl - a2|L|2 n(x) , (2) 

where a is the lattice parameter. In the continuum limit 
(a — ^ 0), this becomes 

N(x) = aL(x) + n(x) + Oia^) . (3) 

Notice that we always have |n(x)p = 1. 

In terms of these and of the continuum fermion field 
?Aq(x) corresponding to Cia we can express the partition 
function as the functional integral 



X exp 



, (4) 
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where the continuum hamiltonian density reads 

U = V^(^(7•V-^)v + Ps|Vn|2+x±52|L|2 

+S'L • [Jk s + i (n X 3rn)] 
+ {-l)\^\SJ^n-s, 



(5) 



with ps = JS^ as the spin stiffness, x± = 4J the trans- 
verse susceptibility and s the itinerant spin operator, given 

by 

S = V'l {^)s-y V'7 ■ (6) 

In expression ([5]), the first term is the continuum elec- 
tron kinetic hamiltonian density, derived from the tight- 
binding energy assuming the system has a Dirac-like dis- 
persion relation near the Fermi points. As the system 
is doped, charge carriers are added or removed from the 
conduction band. Their total number is controlled by a 
chemical potential /i, which has, therefore, been included 
in the previous equation. 

Integrating over L in ([4]) we obtain the resulting effective 
lagrangian density 



+JKl(-l)l"l5n 



Ps 



Vnl' 



1 



\dtn\ 



i , ^ , , Jk 
— (n X OTn) + s 



(7) 



where (7°)^ = 1, 7^7* = cr' and c — y/J>7x± is the spin- 
wave velocity. 

It will be convenient to use the CP^ (Schwinger Boson) 
formulation of the 0(3) NLSM, in which the AF fluctua- 
tion field is written as 



(8) 



in terms of the two complex fields Zq,, a = 1,2, satisfying 
the constraint |zip 4- |z2p = 1. In the CP"'^ language the 
effective lagrangian density ([7]) is rewritten as 



+Jk V'^ 



(-1)1^15 a - n + — (n x d^n) 

X± 
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s • s , 



(9) 



where the components of n are given by ([5]) and 13^ = 

A Gauge Coupling Replaces the Magnetic Inter- 
action . — We now perform a canonical transformation 
[3T| on the electron field, namely 



(10) 



where the unitary matrix U is written in terms of the Za- 
fields as 



U 



zi 

Z2 



(11) 



This matrix has the following property 

U^a-nU = a' (12) 

and therefore the first term in the second line in ^ can 
be expressed, up to a sign, as the density difference of 
electrons with opposite spins: 



1 



,|x| 



JkS (V'tV't - "^Iv-i 



(13) 



Assuming a uniform density of electrons, we conclude 
that this term will vanish upon integration in x because 
of the rapidly oscillating pre- factor. 

Now, since U represents a local operation, it follows 
that, under the transformation pUj) . the electron kinetic 
term generates the additional interaction 



From pT|) . we obtain 



(14) 







-Z2 dnZi + Zi df,Z2 







,(15) 



where we used the fact that = —iz* df^Zi , which follows 
from dni). 

Now, consider the polar representations of the fields Zq,, 



_ Pa iOc — 1 9 

V2 



(16) 



Integration over the di fields eliminates rapidly oscillat- 
ing phase dependent terms, such as the second term in 
P^ . We show in the Appendix that the same happens to 
the second term in the second line in The effective 
lagrangian density, therefore, becomes 



off 



2XJ 



(17) 



Observe that the magnetic interaction between the itin- 
erant electrons and the localized spins manifests ulti- 
mately as a gauge coupling between the electrons and the 
Schwinger boson fields, mediated by the CP^ vector field 
Afj_, which becomes a gauge field. Indeed, (jl7l) is invariant 
under the gauge transformation 



9^ ^ 9, + A, 



A,. 



A^-d^A. 



(18) 



Our aim is to obtain the net effective interaction among 
the conduction electrons, associated to the fermion fields 
Tp. For this purpose we are going to functionally integrate 
over the CP^ fields, in order to derive the resulting ef- 
fective interaction. Before doing that, however we shall 
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express the effective lagrangian in an explicitly gauge in- 
variant way. 

We first introduce the gauge invariant phase-fields |33| 



(19) 
!/□ is the 



which are clearly invariant under (|18p . In (|T! 
Green function of the □ = 9^9^ operator. 

We can now re-write the effective lagrangian density 
(|17p in an explicitly gauge invariant form given by [3 3) 



off 



— ,/,t 



4ps 



-F„, 



72 



-a 



i=l 



s • s . 



(20) 



where we have used the constant pi approximation for i = 
1, 2 and F^^ = dfj.A^ - d^A^. 

The Effective Electron Interaction. Let us now 

perform the functional integration over the bosonic fields 
p, X and Afj_. This will ultimately generate the final effec- 
tive interaction among the conduction electrons. 

We shall adopt the constant pi {i — 1,2) approxima- 
tion for performing the functional integration over the 
Schwinger boson fields, namely, piS and Xi's. This approx- 
imation usually reproduces the physical situation found 
in many materials. It implies that integration over the 
Zi fields would just yield a trivial multiplicative constant 
in the partition function. The nontrivial interaction ef- 
fect comes from integration over the gauge field A^. This 
can be easily performed given the quadratic dependence 
of (Uni) in this field. 

The resulting effective lagrangian density for the con- 
duction electrons is 



s • s 



(21) 



Explicitly writing the components of the Dirac field we 
may, after some algebra, express the effective interaction 
term above as 



4p, 



''2X±' 



s • s . 



(22) 



The first term above is a superconducting, s-wave BCS- 
type interaction. It would lead to a superconducting phase 



with an order parameter A — ('0it'/'2i "t" ''/'2tV'i4,) ■ The sec- 
ond term is a Nambu-Jona-Lasinio-type interaction, that 
may produce an insulating charge-gapped phase showing 
an excitonic condensate. The gap parameter would be 
M — {tptp). Finally, the second and fourth terms in (|22p 
are, respectively, Ising and Heisenberg-like, magnetic in- 
teractions. These may potentially lead to a magnetically 
ordered phase with the nonzero order parameter being the 
average magnetization m = (s). 

A model containing the interactions of the first and 
third terms in (|22l) has been investigated in a mean-field 
approximation [25]. We have shown that, even if the ex- 
citonic interaction strength is larger than the supercon- 
ducting interaction, as the chemical potential increases, 
superconductivity eventually suppresses the excitonic or- 
der parameter, which means that the system goes from an 
insulating state to a superconducting one as charge car- 
rriers are added to the system. 

A model with Dirac fermions subject to the interaction 
described by the first term in (12^ was studied in [T^[T^ . 



Conclusions. — We have provided a concrete exam- 
ple of a model containing both localized spins and itiner- 
ant electrons, where different competing interactions are 
generated out of the original purely magnetic interactions. 
The localized spins present an AF Heisenberg interaction 
on a square lattice, the itinerant electrons are Dirac-like 
and the mutual localized-itinerant interaction is magnetic, 
Kondo-like. 

Our results yield an unified, controllable picture for the 
the common magnetic origin of three competing types of 
order in a strongly correlated system: superconducting, 
magnetic and charge ordering. 

Appendix. — Let us show here that the second term 
in the second line in ([9]) only contains rapidly oscillating 
phase dependent terms, which are, therefore, eliminated 
through functional integration over the phase fields. 

We may write this term as 



{ij^'aip) ■ (n X drii) = ip'lipfjz*^z^dr {z*xZp) 



^^ik „i „i „k 



(23) 



where we have used 

Now, we have the following identity for Pauli matrices 



+16'^'' [-5^^^^, - Sc^a''^,] . (24) 

Inserting (|24p in we immediately obtain 

s • (n X drti) = ipl^il^aZ^Zf^dr {zlzx) + pdt, (25) 

where pdt stands for "phase dependent terms" . 

The first term on the r.h.s. above vanishes because 
z\z\ = 1 and, therefore, only phase dependent terms are 
left, as we have asserted. 
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